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an optimum excitation signal based on the model (207). An application to ultrasonic imaging for 
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Optimization of Excitation Waveform for Nonlinear 
Transmit-Receive Systems 



Background of Invention 
[0001] Field of Invention 

[0002] This invention relates to transmit-receive systems and, more particularly, to 
optimizing the excitation waveform that is used in such systems to maximize the 
output for a given constraint on the input. 

[0003] Description of Related Art 

[0004] Systems employing a transmitter and receiver are widely used. 
Applications include telemetry, imaging and communications. A specific example of 
an imaging application is ultrasonic transmission-mode tomography used in breast 
mammography. 

[0005] A signal of a selected waveform is typically delivered to the transmitter, 
emitted by the transmitter, propagated though a medium, and received at the 
receiver. Information of interest is then extracted. 

[0006] To maximize the signal-to-noise ratio, it is often desirable to maximize the 
output at the receiver based on a given constraint on the input to the transmitter. 

[0007] One approach has been to optimize the design of the individual hardware 
components, such as the transducers and electronics. For example, the thickness 
and impedance of the various matching layers in the transducers have been 
adjusted. Appropriate tuning circuitry has also been introduced to maximize the 
effect of the excitation waveform on the transmitted transducer. 

[0008] These types of tuning methods can be problematic. The impedance of the 
electronic circuits associated with the transmitter and receiver may be different, • 
requiring different tuning conditions during transmission and reception. Precise 
knowledge of the dynamic characteristics of the transducers may also be needed in 
conjunction with the associated electronics. 

[0009] Rather than physically tuning system components, transmit-receive 
systems have been modeled as linear systems. Based on this modeling, the center 
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frequency, shape, duration and amplitude of the excitation waveform has been 
adjusted to maximize the peak received value predicted by these models. 

[0010] Unfortunately, there are often inherent nonlinearities in transmit-receive 
systems, including in ultrasonic transmit-receive systems, that can cause the 
results of this linear modeling approach to be less than optimal. In turn, the greatest 
possible signal-to-noise ratio may not be achieved. 

Summary of Invention 

[0011] A process for optimizing the excitation waveform that is delivered to an 
ultrasonic transmitter that, together with an ultrasonic receiver, form part of a 
nonlinear ultrasonic transmission and reception system. A transmission test signal 
may be delivered to the ultrasonic transmitter. A received test signal may be 
generated from the ultrasonic receiver that is a nonlinear function of the 
transmission test signal. A nonlinear model of the nonlinear function may be 
developed from the transmission test signal and the received test signal. An 
optimum excitation signal for the ultrasonic transmitter may be determined that 
substantially maximizes the signal generated by the ultrasonic receiver based on 
the model and based on a specified constraint on the excitation signal. 

[0012] The kernel functions of the nonlinear function may be determined as part 
of the process of developing the model. An algorithm may be used in determining 
the kernel functions. Principal dynamic modes of the nonlinear function may also be 
determined as part of this process and may be based on the kernel functions. 

[0013] The time inversion of one or more of the principal dynamic modes may be 
determined as part of determining the optimum excitation signal. 

[0014] One or more of the kernel functions may be excluded when calculating the 
time inversion. 

[0015] A Laguerre-Volterra network may be used in developing the nonlinear 
model. Parameters of the Laguerre-Volterra network may be adjusted to minimize 
the mean-squared error between the signal predicted by the network and the 
received test signal. The adjustment may be an iterative process. 

[0016] The nonlinear model may include a linear filter followed by a static 
nonlinearity. 
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[0017] The specified constraint on the excitation signal may include a constraint 
on the amplitude or power of the excitation signal. 

[0018] The amplitude or power of the signal generated by the ultrasonic receiver 
may be maximized as part of determining the optimum excitation signal. 

[001 9] The transmission test signal may be a wideband signal. The wideband 
signal may cover the bandwidth over which the ultrasonic transmitter is configured 
to operate. The wideband signal may be white noise or a chirp. 

[0020] The transmission test signal may cover the dynamic range over which the 
ultrasonic transmitter is configured to operate. 

[0021] The ultrasonic transmitter may be excited with the optimum excitation 
signal. A breast may be placed between the ultrasonic transmitter and the 
ultrasonic receiver during this process. The signal received by the ultrasonic 
receiver in response to the optimum excitation signal may be analyzed to create an 
image of the breast. 

[0022] An ultrasonic imaging system for generating an image of tissue may 
include an ultrasonic transmitter that converts an excitation signal into an ultrasonic 
signal; an ultrasonic receiver positioned to receive the ultrasonic signal transmitted 
by the ultrasonic transmitter and that generates a received signal that is a nonlinear 
function of the excitation signal; an excitation signal generator in communication 
with the ultrasonic transmitter that generates an excitation signal that substantially 
maximizes the signal generated by the ultrasonic receiver based on a specified 
constraint on the excitation signal; and a processing system in communication with 
the ultrasonic receiver for processing the signal generated by the ultrasonic 
receiver into an image of tissue disposed between the ultrasonic transmitter and the 
ultrasonic receiver. 

[0023] The signal generated by the excitation signal generator may be derived 
from a nonlinear model of the nonlinear function. The nonlinear model may be 
developed from a comparison of a transmitted test signal transmitted by the 
ultrasonic transmitter and a received test signal generated by the ultrasonic 
receiver. 
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[0024] These processes and system may also be used in transmit-receive 
systems that operate outside of the ultrasonic range. 

[0025] These as well as still further features, benefits and objects will now 
become clear upon an examination of the following detailed description of 
illustrative embodiments and the attached drawings. 

Brief Description of Drawings 

[0026] Fig. 1 is a block diagram of a transmit-receive system. 

[0027] Fig. 2 is a block diagram of a process for determining the optimum 
excitation signal for a transmitter in a non-linear transmit-receive system, 

[0028] Fig. 3 is a diagram of a Laguerre Volterra Network. 

[0029] Fig. 4 is a diagram of parallel cascades of linear filters and their associated 
nonlinearities. 

[0030] Fig. 5 is non-linear model of a transmit-receive system using a linear filter 
and a static nonlinearity. 

J0031] Figs. 6(a) - 6(d) are traces of measured and predicted signals in a non- 
linear transmit-receive system. 

[0032] Fig. 7 is a block diagram of a non-linear model of a non-linear system 
using principal dynamic modes. 

[0033] Figs. 8(a) - (d) are traces of optimal and pulse excitation signals and their 
corresponding received signals. 

Detailed Description of Illustrative Embodiments 
[0034] Fig. 1 is a block diagram of a transmit-receive system. 

[0035] As shown in Fig. 1, an excitation signal generator 101 may generate an 
excitation signal. This signal may be a steady-state signal, a series of pulses, or 
any other type of signal. 

[0036] The signal generated by the excitation signal generator 101 may be 
delivered to a transmitter 103. The transmitter 103 may be an ultrasonic transmitter 
or any other type of transmitter. The transmitter may be a single transmitter or an 
array of transmitters. In the case of an array of transmitters, separate signals may 

-4- 

5/3/05, EAST Version: 2.0.1.4 



WO 2004/010590 



PCT/US2003/022531 



be delivered to each element in the array from the excitation signal generator 101 . 
The separate signals may be the same or different. 

[0037] The signal transmitted by the transmitter 103 may travel through a 
propagation medium, such a breast suspended within acoustical coupling fluid, and 
then arrive at a receiver 105. The receiver 105 may be an ultrasonic receiver or any 
other type of receiver. As with the transmitter 103, the receiver 105 may be a single 
receiver or an array of receivers. 

[0038] The signal generated by the ultrasonic receiver 105 (or signals in the case 
of an array) may be directed to a processing system 107. The processing system 
may analyze the received signal, possibly along with the excitation signal, to extract 
information that is of interest. In one embodiment, the information of interest may . 
be a three-dimensional image of a breast interposed between the transmitter 103 
and the receiver 105. More details concerning transmission mode tomographic 
imaging for breast mammography are set forth in U.S. Patent Application Serial No. 
10/117,763, filed April 5, 2002, entitled "High Resolution 3D Ultrasonic 
Transmission Imaging," U.S. Patent Application Serial No. 10/117,468, filed April 5, 
2002, entitled "Nonlinear Processing for Mitigation of Diffraction Effects," U.S. 
Provisional Application Serial No. 60/396,516, filed July 16, 2002, entitled "Support 
Bra for Ultrasonic Breast Scanner, " U.S. Patent Application Serial No. 10/382,235, 
filed March 5, 2003, entitled "Multi-Mode Processing for Ultrasonic Imaging," U.S. 
Provisional Application Serial No. 60/478,981, filed June 17, 2003, entitled 
"Targeted Ultrasonic Hyperthermia With A High-Resolution 3-D Fixed Coordinate 
System," and U.S. Patent Application filed June 27, 2003, entitled "Three- 
Dimensional Ultrasonic Mammography Scanning," attorney docket number 64693- 
070. The entire content of all of these applications is incorporated herein by 
reference. 

[0039] Fig. 2 is a block diagram of a process for determining the optimum 
excitation signal for a transmitter in a nonlinear transmit-receive system. 

[0040] The transmitter 103 may be excited with a test signal, as reflected by an 
Excite Transmitter With Test Signal step 201 . The test signal may be a wideband 
signal. It may be bandwidth-limited white noise, such as white noise with a 
bandwidth of 4 - 24 MHz. It may be a chirp. The wideband test signal may be 
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selected to cover the bandwidth and dynamic range over which the transmitter 103, 
the receiver 105 and their associated electronics are configured to operate during 
normal use. The test signal may be any other form of signal. 

[0041] The test signal may travel through the anticipated propagation medium, 
such as acoustic coupling fluid containing a breast dangling therein. The test signal 
may then reach the receiver 105. The receiver 105 may then generate a received 
test signal, as reflected by a Receive Test Signal from Receiver step 203. 

[0042] The excitation signal and the signal generated by the receiver may be 
stored for later analysis. If desired, the signals may be digitized either before or 
after their storage. A time-bandwidth product of approximately 1 ,000 may be used 
for this purpose which, in the case of an excitation signal having a bandwidth of 
approximately 4-24 Megahertz, may be approximately 100 MHz. 

[0043] The signal generated by the receiver 105 may be a nonlinear function of 
the excitation signal generated by the excitation signal generator 101. The 
nonlinearity may be caused by nonlinearity in the transmitter 103, the receiver 105, 
the excitation signal generator 101, the propagation medium, and/or by other areas 
of the system. 

[0044] A nonlinear model of the dynamic transformation between the excitation 
signal and the signal generated by the receiver may be developed, as reflected by 
a Develop Nonlinear Model step 205. This may be a mathematical model that is 
derived using any type of nonlinear system identification method. 

[0045] Based on the nonlinear model that is developed, an optimum excitation 
signal may be determined that optimizes the signal generated by the receiver 105 
based on a constraint on the excitation signal that is delivered to the transmitter 
103. Any technique for determining the optimum excitation signal based on the 
nonlinear model may be used. 

[0046] The constraint on the excitation signal may include a constraint on its peak 
amplitude, power and/or other characteristics. 

[0047] The excitation signal may be optimized to maximize the amplitude, power 
and/or other characteristics in the signal generated by the receiver 105. 
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[0048] The optimum excitation signal may then be applied to the transmitter 103, 
as reflected by an Apply Optimum Excitation Signal step 209. To this end, the 
excitation signal generator 101 may be configured to generate the optimum 
excitation signal during the normal operation of the transmit-receive system. 

[0049] Fig. 3 is a diagram of a Laguerre-Volterra network (LVN). This is an 
example of a nonlinear system identification method that may be used as part of 
the process of developing a nonlinear model of the transmit-receive system. 

[0050] In Fig. 3, x(n) represents the excitation signal generated by the excitation 
signal generator 101 and y(n) represents a prediction of the signal that is generated 
by the receiver 105 as a result. 

[0051] The LVN method may combine expansions on the discrete-time Laguerre 
basis with feed-forward artificial neural networks (FANN) using polynomial 
activation functions. It may follow the standard architecture of a single-layer, fully 
connected FANN with three distinctive features: 

[0052] (1 ) The Laguerre-filter bank bj (/=0, 1 , . . . , L) that preprocesses the 

input, x(n); 

[0053] (2) The polynomial activation functions, f k , (fc=1 ,2, ... v K) in the 

hidden units; and 

[0054] (3) The non-weighted summative output unit. 

[0055] The architectural parameters of the LVN may consist of the number of 
Laguerre filters (LFs), hidden units (HUs) and the degree of the polynomial 
activation functions. The selection of these structural parameter values can be 
performed by successive trials in ascending order, moving from lower to higher 
numbers, using the normalized mean-squared error (NMSE) of the output 
prediction achieved by the model as a criterion. 

[0056] The determination of the optimal number of LFs for a specific system can 
be assisted by observing the relative magnitude of the Laguerre expansion 
coefficients. Likewise, the determination of the optimal number of HUs can be 
assisted by observing the relative magnitudes of the polynomial coefficients, when 
the inbound weights for each HU are normalized to a unity sum of squares. 
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[0057] Once the LVN architecture is set, the training of the LVN may be 
performed with input-output experimental data to estimate the values of the 
unknown network parameters, including the weights Wkj, the polynomial coefficients 
c q ,k, the Laguerre parameter a and the offset y 0 . The network training may involve 
minimizing the mean-square error between the measured output and the LVN 
predicted output y(n) for the corresponding input. For effective training, the input- 
output data may be broadband and cover the entire bandwidth and dynamic range 
of the system. 

[0058] The output of the y-th discrete-time Laguerre filter with impulse response 
function bj(m) may be given by the discrete convolution: 

[0059] vj(n) = %bj(m)x(n-m) (1) 

[0060] where 

[0061] bJm) = a™''(l-a)'''£H/ . J a^'(l-a)' (2) 

i=o \ l A 1 J 

[0062] and a is the discrete-time Laguerre parameter. 

[0063] The input of the k-th hidden unit may be a weighted sum of the Laguerre 
filter bank outputs: 

[0064] u k (n) = Jv/^ (k=l,2,..K) (3) 

[0065] The corresponding output of the Mi hidden unit may be given by the 
polynomial activation function: 

t 

q-0 



[0066] z k (n)=2Lc qM u\(n) (k = 1.2,...K) (4) 



[0067] The network output may be given by a non-weighted summation of the 
hidden unit outputs including a trainable offset value, yo: 

[0068] y(n) = f i Z k (n) + y 0 (5) 
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[0069] The training of the network parameters may be performed iteratively using 
a gradient descent method for a quadratic cost function defined by the square of 
the output prediction error. 

[0070] As shown in Fig. 4, the LVN architecture shown in Fig. 3 may generally be 
expressed as K parallel cascades of linear filters ... L k and their associated static 
nonlinearities A/* . . . N k . Fig. 4 shows the equivalent version of the LVN where the 
input x(n) goes through a set of K linear filters (determined by the weights w kJ ) and 
their associated nonlinearities (determined by the polynomial coefficients c qtk ). The 
output may be formed by summing all of the outputs of the static nonlinearities. 

[0071] More details concerning the use of a Laguerre-Volterra network as a 
nonlinear system identification method are set forth in V.Z. Marmarelis, "Modeling 
Methodology for Nonlinear Physiological Systems," Annals of Biomedical 
Engineering, Vol. 25, pp. 239-251 (1997), the entire content of which is 
incorporated herein by reference. 

[0072] Once the network parameters are found, the corresponding model of 
kernel functions may be constructed in accordance with well known techniques. 
The obtained mathematical model may describe the nonlinear dynamic 
transformation of the input excitation signal into the output signal generated by the 
receiver in the transmit-receive system. 

[0073] In one experiment, the test signal was a wideband chirp covering a 
bandwidth from approximately 4 to 20 MHz over a dynamic range of ± 50 volts. The 
received signal was recorded over 5,000 samples at 20 nsec sampling intervals. 

[0074] The excitation waveform was generated from an arbitrary waveform 
generator, such as an HP 33250Afrom Hewlett-Packard in Palo Alto, California, 
with a pre-designed wideband waveform covering this bandwidth and dynamic 
range. The system output was measured at the receiving transducer and digitized 
using an oscilloscope, such as the TDS5054 Digital Phosphor Oscilloscope from 
Tektronix in Beaverton, Oregon. 

[0075] The digitized transmitted and received signals were analyzed using the 
LVN nonlinear system identification method described above to obtain a nonlinear 
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model of the dynamic transformation of this transmit-receive system for any 
excitation waveform into the corresponding received waveform. 

[0076] Application of the LVN modeling methodology to the collected data yields 
the nonlinear model shown in Fig. 5. This is comprised of the excitation waveform 
501 being delivered to a linear filter 503, followed by a cubic static nonlinearity 505 
and a received signal 507. 

[0077] With the LVN structural parameters of 17 Laguerre functions, 2 hidden 
units, and a third-degree polynomial activation function, the normalized mean 
squared-error of the model prediction with respect to the measured output data was 
5.9% after 100 iterations. In all trials with different structural parameters, only one 
hidden unit was found to be active. This allowed the model to be simplified to the 
single cascade shown in Fig. 5. 

[0078] Fig. 6 shows the measured receiver output and the LVN model predictions 
along with their power spectra. The accuracy of the LVN model prediction is evident 
in both the time and the frequency domains. 

[0079] The optimal excitation waveform may be derived from this model on the 
basis of the matched-filter principle, that is, by determining the time-reversed 
impulse response of the linear filter of the cascade model, scaled according to an 
input power constraint. The optimal excitation pulse may match the frequency 
response of the linear filter. The static nonlinearity may determine the final gain of 
the output and indicate that doubling of the input more than doubles the output 
(supralinear relation). 

[0080] Fig. 7 is a block diagram of a nonlinear model of a nonlinear system using 
principal dynamic modes. The Pf . . .P k blocks represent the principal dynamic 
modes (PDM) of the system based on the kernel functions. The ui(t) . . . u k (t) 
represent a convolution of the input xft) with each of the principal dynamic modes. 
The f block represents the static point-to-point mapping. 

[0081] Mathematically, y(t) in Fig. 7 may be expressed as: 

[0082] y(0 = fk(t),.~,u k (t)} (6) 

[0083] The input signal x(t) may be expressed as: 
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[0084] *(0 = 2>/ P /K> (?) 

[0085] where Ai are scalers. 

[0086] By applying the input signal, uj(t) may be expressed as: 

[0087] Uj (t) = ^ A, j/? y (j)p { (r-tytr (8) 

M 0 

[0088] =2]W0 (9) 

[0089] and y(t) may be expressed as: 

[0090] ^0 = /{V-A} (10) 

[0091 ] where Zy (JJ are known. 

[0092] yft) may then be optimized over the parameters {A ; ,...,A A }subject to a 
constraint: 

[0093] £<a; X k ) (11) 

[0094] For example, the values (Ai Ak) may be determined for the maximum 

value of y(t) over t is greatest subject to the constraint: 

[0095] %(K) 2 = constant (12) 

[0096] Other optimization criteria, such as maximizing the power of y(t) and/or 
constraints such as 0 < Ai < 1 , can also be applied. 

[0097] As indicated, these are examples of approaches for determining an 
optimum excitation signal based on an input constraint. 

[0098] Figs. 8(a) - (d) are traces of optimal and pulse excitation signals and their 
corresponding received signals. Fig. 8(a) is a trace of the optimal excitation signal 
that was determined based on the principals discussed above; Fig. 8(b) is a trace 
of the signal generated by the receiver of this optimal excitation signal; Fig. 8(c) is a 
trace of a non-optimal pulse that was commonly used to excite a transmitter; and 
Fig. 8(d) is a trace of the signal generated by the receiver of this non-optimal pulse 
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signal. As can be seen by a comparison of the traces in Figs. 8(b) and (d), use of 
the optimal excitation signal increased the amplitude of the received signal by 
about 35 times, without any significant increase in the amplitude of the input signal. 

[0099] Although a test signal has been described as a mechanism for obtaining 
the nonlinear model of the system, other approaches could be used instead. In 
short, what has thus-far been described are merely examples of the components, 
steps, features and benefits of the invention. The invention is not limited to these, 
but solely to the subject matter delineated by the following claims and their 
equivalents. 

[0100] This application is based upon and claims priority to U.S. Provisional 
Application Serial No. 60/397,378, entitled "Optimization of the Excitation Waveform 
for Transmit-Receive Systems," filed July 19, 2002, the entire content of which is 
incorporated herein by reference. 
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We claim: 

1 . A process for optimizing the excitation waveform that is delivered to 
an ultrasonic transmitter that, together with an ultrasonic receiver, form part of a 
nonlinear ultrasonic transmission and reception system comprising: 

delivering a transmission test signal to the ultrasonic transmitter; 

generating a received test signal from the ultrasonic receiver that is a 
nonlinear function of the transmission test signal; 

developing a nonlinear model of the nonlinear function from the 
transmission test signal and the received test signal; and 

determining an optimum excitation signal for the ultrasonic transmitter 
that substantially maximizes the signal generated by the ultrasonic receiver 
based on the model and based on a specified constraint on the excitation 
signal. 

2. The process of Claim 1 wherein developing the nonlinear model 
includes determining kernel functions of the nonlinear function. 

3. The process of Claim 2 wherein an algorithm is used in determining 
the kernel functions. 

4. The process of Claim 2 wherein developing the nonlinear model 
includes determining principal dynamic modes of the nonlinear function based on 
the kernel functions. 

5. The process of Claim 4 wherein determining an optimum excitation 
signal includes calculating the time inversion of one or more of the principal 
dynamic modes. 

6. The process of Claim 5 wherein one or more of the kernel functions 
are excluded when calculating the time inversion. 

7. The process of Claim 1 wherein developing the nonlinear model 
includes determining principal dynamic modes of the nonlinear function. 

8. The process of Claim 1 wherein a Laguerre-Volterra network is used 
in developing the nonlinear model. 
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9. The process of Claim 8 wherein parameters of the Laguerre-Volterra 
network are adjusted to minimize the mean-squared error between the signal 
predicted by the network and the received test signal. 

10. The process of Claim 9 wherein the adjustment is an iterative 

process. 

1 1 . The process of Claim 1 wherein the nonlinear model includes a linear 
filter followed by a static nonlinearity. 

12. The process of Claim 1 wherein the specified constraint on the 
excitation signal includes a constraint on the amplitude of the excitation signal. 

13. The process of Claim 1 wherein the specified constraint on the 
excitation signal includes a constraint on the power of the excitation signal 

14. The process of Claim 1 wherein determining an optimum excitation 
signal maximizes the amplitude of the signal generated by the ultrasonic receiver. 

15. The process of Claim 1 wherein determining an optimum excitation 
maximizes the power of the signal generated by the ultrasonic receiver. 

16. The process of Claim 1 wherein the transmission test signal is a 
wideband signal. 

17. The process of Claim 16 wherein the wideband signal covers the 
bandwidth over which the transmission and reception system is configured to 
operate. 

18. The process of Claim 16 wherein the wideband signal is white noise. 

19. The process of Claim 16 wherein the wideband signal is a chirp. 

20. The process of Claim 1 wherein the transmission test signal covers a 
dynamic range. 

21 . The process of Claim 20 wherein the dynamic range is the dynamic 
range over which the transmission and reception system is configured to operate. 

22. The process of Claim 1 further comprising exciting the ultrasonic 
transmitter with the optimum excitation signal. 
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23. The process of Claim 22 wherein a breast is placed between the 
ultrasonic transmitter and the ultrasonic receiver while exciting the ultrasonic 
transmitter with the optimum excitation signal. 

24. The process of Claim 23 wherein the signal received by the ultrasonic 
receiver in response to the optimum excitation signal is analyzed to create an 
image of the breast. 

25. An ultrasonic imaging system for generating an image of tissue 
comprising: 

an ultrasonic transmitter that converts an excitation signal into an 
ultrasonic signal; 

an ultrasonic receiver positioned to receive the ultrasonic signal 
transmitted by the ultrasonic transmitter and that generates a received 
signal that is a nonlinear function of the excitation signal; 

an excitation signal generator in communication with the ultrasonic 
transmitter that generates an excitation signal that substantially maximizes 
the signal generated by the ultrasonic receiver based on a specified 
constraint on the excitation signal; and 

a processing system in communication with the ultrasonic receiver for 
processing the signal generated by the ultrasonic receiver into an image of 
tissue disposed between the ultrasonic transmitter and ultrasonic receiver. 

26. The system of Claim 25 wherein the signal generated by the 
excitation signal generator is derived from a nonlinear model of the nonlinear 
function. 

27. The system of Claim 26 wherein the nonlinear model is developed 
from a comparison of a transmitted test signal transmitted by the ultrasonic 
transmitter and a received test signal generated by the ultrasonic receiver. 

28. A process for optimizing the excitation waveform that is delivered to a 
transmitter that, together with a receiver, form part of a nonlinear transmission and 
reception system comprising: 

delivering a transmission test signal to the transmitter; 
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generating a received test signal from the receiver that is a nonlinear 
function of the transmission test signal; 

developing a nonlinear model of the nonlinear function from the 
transmission test signal and the received test signal; and 

determining an optimum excitation signal for the transmitter that 
substantially maximizes the signal generated by the receiver based on the 
model and based on a specified constraint on the excitation signal. 

29. A nonlinear transmission and reception system comprising: 

a transmitter that converts an excitation signal into a transmitted 
signal; 

a receiver positioned to receive the transmitted signal and that 
generates a received signal that is a nonlinear function of the excitation 
signal; and 

an excitation signal generator in communication with the transmitter 
that generates an excitation signal that substantially maximizes the signal 
generated by the ultrasonic receiver based on a specified constraint on the 
excitation signal. 

30. The system of Claim 29 wherein the signal generated by the 
excitation signal generator is derived from a nonlinear model of the nonlinear 
function. 

31. The system of Claim 30 wherein the nonlinear model is developed 
from a comparison of a transmitted test signal transmitted by the transmitter and a 
received test signal generated by the receiver. 
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